SALAMANDRA 48(3) 


ISSN 0036-3375 


157-165 30 October 2012 


Age and size at metamorphosis of Salamandra infraimmaculata larvae 
born in the laboratory and raised under different density regimes 

with food ad libitum 


Michael R. Warburg 


Department of Biology, Technion - Israel Institute of Technology, Haifa 32000, Israel 
e-mail: warburg@tx.technion.ac.il 
Manuscript received: 29 November 2011 


Abstract. Age and size at metamorphosis were studied in two half-sibling Salamandra infraimmaculata larval cohorts 
born and raised in the laboratory under three different density regimes, and fed ad libitum until they metamorphosed. 
There was no significant effect of density on the number of larvae metamorphosing. In the two cohorts studied, a signifi¬ 
cant positive relationship was observed with age at metamorphosis that increased with density, and a significant negative 
relationship between mass and density. In both cohorts, density did not appear to have an effect on either minimal or 
maximal age, mass or length at metamorphosis, nor on the range between maximum and minimum. There was no signifi¬ 
cant difference between the two cohorts in either age or length at metamorphosis. The difference in mass increased with 
density. The evolutionary significance of density effects on size (mass, length) and age at metamorphosis under unlimited 
resource (food) conditions, is by spreading out emergence of post-metamorphs onto land and their subsequent dispersal, 
and by affecting their size as adults and thereby their eventual maturity. 
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Introduction 

Urodeles exhibit plasticity in the timing of, and in achiev¬ 
ing an optimal size at metamorphosis (Wilbur & Collins 
i973> Warburg 2011). The duration of the larval period 
(i.e., their age or the time to metamorphosis), as well as 
their size at metamorphosis depend on their larval growth 
history (discussed in Cohen et al. 2005, 2006). Both the 
timing of, and the size at, metamorphosis are critical for 
survival of the young post-metamorphs and as a conse¬ 
quence, for the recruitment of juveniles into the terrestrial 
population (Warburg 2011). 

This is especially true when conditions are suboptimal 
as are those prevailing on Mt. Carmel, Israel, during the 
breeding season of Salamandra infraimmaculata Mar¬ 
tens, 1885. Winter rains may start in October, eliciting 
breeding in the salamanders. In some years, rainfalls may 
be sufficient to fill the rock pools, which are their main 
breeding sites. However, in about 50% of the time, the pools 
dry out, resulting in the loss of entire larval cohorts (War¬ 
burg 2009a). Nevertheless, when rains eventually resume, 
some larvae may be saved, since they can survive for about 
a week on wet mud (Warburg 1986a,b). 

The question remains: why start breeding early in the 
rainy season when about half of the broods might be lost? 
The advantage in starting larval life as early as possible in 


the season, lies in the advantage these larvae have when 
rainfall is continuous, since they are able to grow more 
rapidly by feeding on newborn larvae (Cohen et al. 2005, 
Warburg 1992, 2009a, 2010). The larvae require at least 
five weeks to undergo metamorphosis (Warburg 2011). If 
breeding is delayed because of delayed rainfall, there is a 
risk that time could be insufficient to allow metamorphosis 
(Warburg 2010, 2011). 

Aquatic fauna is rather poor when the ponds first form 
early in the season. Larvae then depend almost entirely on 
their cannibalistic (sibling predation) trait, preying on late- 
born larval cohorts in order to obtain enough food to facil¬ 
itate their own growth (Cohen et al. 2005). The survival of 
larvae in spite of all these risks will enable their successful 
survival as post-metamorphs. 

In addition, the emergence on land of the young post- 
metamorphs will have to be timed (i.e., actual date) so that 
the soil will still be moist and temperatures sufficiently low 
to make successful dispersal possible (Warburg 2011). 
This depends on the amount of rains falling until then, and 
on the season. Emergence in winter is preferable to spring, 
which is characterized by hot spells (sciroccos). Moreover, 
both the larvae’s age at metamorphosis (measured from 
their birth in the laboratory because it can rarely be re¬ 
corded in nature), as well as the duration of the larval pe¬ 
riod (i.e., the time it takes to metamorphosis), are of great 
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significance. As a rule, larval period is positively related to 
age at metamorphosis. Finally, the size attained by the lar¬ 
vae when metamorphosing is of great consequence for the 
survival of the post-metamorphs and later of the adults. 
However, metamorphosis can be delayed due to slow lar¬ 
val growth rates, resulting in small size at metamorphosis. 
To conclude: successful survival to adulthood depends on 
starting life as a larva as early in the season as made pos¬ 
sible, by rainfall (Warburg 1992, 2010), growing rapidly 
(Cohen et al. 2006), avoiding cannibalistic predation by 
half-siblings (Cohen et al. 2005), metamorphosing as ear¬ 
ly in the season as possible and at an early age (Warburg 
2011), attaining a greater size, and, finally, dispersing as ear¬ 
ly in the season as possible in order to escape dehydration 
of the soil. 

Density (i.e., larval crowding) was shown to have an ef¬ 
fect on urodelan larvae with regard to both growth rate 
and development to metamorphosis (Brunkow & Col¬ 
lins 1996). Ohdachi (1994) has shown that density was 
negatively related to both time and size at metamorphosis. 
High density can cause failure in metamorphosis (Wilbur 
1976). Scott (1990) in Ambystoma opacum (Graven- 
horst, 1807), and Brodman (1996) have shown in both 
Ambystoma maculatum (Shaw, 1802) and A. jeffersonia- 
num (Green, 1827) that both size and larval period are 
negatively affected by density. 

The objective of the present study is to establish to which 
extent three different density regimes affect both age and 
size (mass, length) at metamorphosis of half-sibling sala¬ 
mander larval cohorts born in the laboratory to the same 
female (i.e., a cohort), and fed ad libitum. 

Materials and methods 

The data discussed here are based on a long-term study 
(1974-1998) of a single breeding population of Salaman- 
dra infraimmaculata (Warburg 2006, 2007a). This species 
is a rare’ and endangered species found in the northern, 
mountainous part of Israel (IUCN Red List status “Near 
Threatened”; www.iucnredlist.org). The metapopulation 
inhabiting Mt. Carmel survives in the most southeastern 
portion of the genus’ distribution. 

Gravid females collected in the field delivered their lar¬ 
vae in the laboratory and were thereafter released into the 
wild (Warburg 2007b, 2008, 2009a). This is an ovovivipa- 
rous species, laying eggs that hatch upon contact with water 
(Warburg et al. 1978/79). During the study period, a to¬ 
tal of 74 half-sibling larval cohorts were born in the labora¬ 
tory (Warburg 2009c, 2010), two of which were used for 
the study described in this paper. Thus, the identity of their 
mothers, their dates of birth, and their sizes at birth were 
all known. These cohorts contained 4,085 larvae, all born to 
these freshly collected females. Most of the larvae were re¬ 
leased into the ponds where their mothers had been collect¬ 
ed. Of the remaining larvae, 396 were raised in the laborato¬ 
ry (at room temperature) until they metamorphosed (Co¬ 
hen et al. 2005, 2006). Others were raised for an additional 


two to five years after metamorphosis in order to release 
them into the wild as juveniles (Warburg unpubl. data). 

Larvae were placed in glass ‘finger bowls’ (13.5 cm in di¬ 
ameter) or enamel troughs (15 x 25 cm) filled with stale 
(aged) tap water, 2 cm in depth. They were fed ad libitum 
(i.e., the larvae had available unlimited amounts of food 
until they completed metamorphosis) with live Tubifex 
worms or minced beef liver (no significant differences were 
found in the effects of these two types of food on either age 
or size at metamorphosis). The Tubifex were added after 
cleaning the troughs and were always available to larvae. 
The water was changed every day. 

Larval salamanders usually stay on the bottom most of 
the time, hiding among stones. Only towards the end of 
the metamorphic cycle will they swim to the surface for 
air. Consequently, density is calculated per surface area 
rather than volume. Larvae were placed either individually 
into finger bowls (42.4 cm 2 /larva), or in enamel troughs 
in groups of either five (‘low density’, 28.7 cm 2 /larva) or 10 
larvae (‘high density’, 14.3 cm 2 /larva). Three replicates of 
each constellation were used. Towards the end of their lar¬ 
val period, they were transferred into spacious plastic con¬ 
tainers with both water and soil that enabled them to crawl 
onto land and metamorphose. 

As the sample size was rather small, I used both t-tests 
and regression analysis for statistical analyses with Sigma- 
Plot 9.0. The relationships between age, mass, and length, 
and the three density conditions were tested using regres¬ 
sion analysis. The difference between the two cohorts was 
examined by t-tests. 

Results 

The data for means and ranges in age and dimensions 
(mass, length) of the two half-sibling larval cohorts are giv¬ 
en in Tables 1 and 2. Cohort I consisted of 96 larvae, and co¬ 
hort II of 106 larvae. There was no significant difference be¬ 
tween the two cohorts (Appendix 1). No significant effect of 
densities on the numbers of emerging metamorphs could 
be seen in either cohort (Figs. 1, 2) (Appendix 2). A signifi¬ 
cant positive relationship between increasing age and den¬ 
sity was noted in cohort I (Figs. 3A-C) (Appendix 3). Like¬ 
wise in cohort II, age at metamorphosis related significantly 
to densities as did mass and length (Fig. 4A-C). In both 
cohorts, no significant effect of density could be demon¬ 
strated as affecting either mass or length with age (Figs. 5A- 
C, 6 A-C) (Appendix 4). Density did not appear to have 
an effect on either minimal or maximal age except at high 
densities in cohort II. In both cohorts, minimal and maxi¬ 
mal mass and length did not show any significant difference 
(Fig. 7 A-C). The range between minimal and maximal age 
at metamorphosis showed an insignificant positive rela¬ 
tionship with increasing density in both cohorts (Fig. 8 A), 
whereas both mass and length at metamorphosis decreased 
with density. Finally, there was no significant difference in 
either age, mass or length between metamorphosing larvae 
belonging to the two cohorts (Fig. 9 A-C) (Appendix 1). 
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Table 1. Age, mass, and length at metamorphosis of specimens in cohort I, fed ad libitum and raised under different density regimes 
(single: larvae kept individually in glass jars (42.4 cm 2 /larva); low = larvae kept at a density of 28.7 cm 2 /larva; high = larvae kept at a 
density of 14.3 cm 2 /larva). Means are followed by standard deviation, the range in parentheses and the sample size (N). 



single 

low 

high 

age (days) 

44.1 ± 2.04 (41-47); N = 6 

45.45 ± 2.56 (40-51); N = 6 

46 ± 3.28 (41-76); N = 6 

mass (g) 

1.63 ± 0.14 (1.46-1.6); N = 16 

1.36 ± 0.22 (1.03-1.79); N = 16 

1.29 ± 0.23 (1.07-1.55); N = 16 

length (mm) 

66.3 ± 3.76 (64-70); N = 6 

60.1 ± 1.87 (57-64); N = 16 

60.2 ± 3.16 (56-70); N = 10 


Table 2. Age, mass, and length at metamorphosis of specimens in cohort II, fed ad libitum and raised under different density regimes. 
For abbreviations see Table 1. 



single 

low 

high 

age (days) 

44.4 ± 2.07 (61-70); N = 10 

44.6 ± 4.28 (40-51); N = 9 

47.9 ± 7.91 (41-76); N = 16 

mass (g) 

1.69 ± 0.09 (1.52-1.83); N = 10 

1.51 ± 0.23 (1.24-1.92); N = 9 

1.44 ± 0.21 (1.08-1.97); N = 16 

length (mm) 

65.6 ± 2.95 (61-70); N= 11 

61.8 ± 3.68 (57-68); N = 9 

59.2 ± 3.54 (54-65); N = 16 




Figure 1. Density effects on numbers of metamorphs in cohort I. 
(A = single, B = low, C = high). See Appendix 2. 
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Figure 2. Density effects on numbers of metamorphs in cohort 
II. (A = single, B = low, C = high). See Appendix 2. 
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Age 

y = 1.1 x +43.083 453 





Figure 3. Average (A) age, (B) mass, and (C) length at metamor¬ 
phosis in cohort I when raised under three different densities. 
See Appendix 3. 

Discussion 

The main problem with studies on aquatic amphibian lar¬ 
vae concerns their origin and age, both of which are gener¬ 
ally unknown. Thus, in some studies on density effects on 
metamorphosis in urodeles, neither age nor larval period 
could be known with accuracy since the larvae studied had 
been collected in the field (Chazal et al. 1994). In others, 
individual egg clutches in nests attended by females were 
collected in the field, separated, and hatched in the labora¬ 
tory. Consequently, although the larval period was known 
(since the eggs hatched in the laboratory) in such cases, 
it was nevertheless not identical with age (from the actual 
point of time when the eggs had been laid) (Kusano 1981, 
Collins & Cheek 1983, Fauth, Resetarits & Wilbur 
1990, Doody 1996, Nishihara 1996, Walls 1998). More¬ 
over, in neither scenario can it be ascertained whether the 
larvae were indeed half-siblings, because it cannot be ruled 
out that either the larvae or the batches of eggs indeed orig¬ 
inated from more than one female. In some other studies, 


Age 





Figure 4. Average (A) age, (B) mass, and (C) length at metamor¬ 
phosis in cohort II when raised under three different densities. 
See Appendix 3. 

both age and larval period were known since individu¬ 
al clutches of half-siblings were either collected (Doody 
1996), or adult females laid eggs that hatched in the labora¬ 
tory (Ziemba & Collins 1999, Cohen et al. 2005, 2006). 

Some of the variability found in the period to metamor¬ 
phosis and size of post-metamorphs can be attributed to 
these facts. Thus, it is not surprising that Wilbur (1976) 
noticed a high variance in both time and size at metamor¬ 
phosis within the same pond. Variation in duration of lar¬ 
val period and body size at metamorphosis was attribut¬ 
ed to time of oviposition and density of conspecific larvae. 
Such individual variation can affect population dynamics 
(Brunkow & Collins 1996). 

Most studies attempted to show the effects of density 
(i.e., crowding or sharing the space with one or more lar¬ 
vae), and resource (i.e., food resources) on time required 
until metamorphosis and size of post-metamorphs. Few 
studies tried to isolate and analyse these factors separate¬ 
ly: the effect of density when food is available at all times, 
or the effect of food when larvae are raised singly without 
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Figure 5. Age, mass, and length at metamorphosis in cohort I 
(A: single, B: low density, C: high density). See Appendix 4. 


any potential density effect (Brodman 1999, Alcoben- 
das et al. 2004, Warburg 2009b, 2010). In the last men¬ 
tioned study it was concluded that metamorphosis appears 
to be affected by food resources when density could not 
have been a factor since these larvae were raised individu¬ 
ally. Table 3 summarizes density conditions in six different 
studies of larvae belonging to four different urodelan spe¬ 
cies that were raised on an unlimited food supply. The con¬ 
ditions vary greatly between different experiments. Thus, 
low-density conditions ranged from 16.5 cm 2 per larva of 
Hynobius retardatus (Dunn, 1923) (see Wakahara 1995) 
to 502.7 cm 2 per larva of Ambystoma tigrinum nebulosum 
Hallowell, 1853 (Collins & Cheek 1983): a 30.5-fold dif¬ 
ference. Likewise, high-density conditions ranged from 
7.5 cm 2 /larva of Salamandra salamandra (Linnaeus, 1758) 
(Degani 1993) to 42.7 cm 2 /larva of Ambystoma texanum 
(Matthes, 1855) (petranka 1984): a 5.7-fold difference. 
In the present study, S. infraimmaculata larvae were raised 
at a low density of 28.7 cm 2 /larva or a high density of 14.3 
cm 2 /larva. Since experimental conditions vary so much 
between different studies on the significance of density 
in urodelan larval growth and metamorphosis, it is diffi- 
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Figure 6. Age, mass, and length at metamorphosis in cohort II 
(A = single, B = low density, C = high density). See Appendix 4. 


cult to compare them. Thus, the high density of 42.7 cm 2 
in A. texanum larvae (Petranka 1984), is close to a mod¬ 
erate density of 43.7 cm 2 in S. salamandra larvae (Dega¬ 
ni 1993). Raising larvae of A. maculatum at high densities 
significantly affected the number of metamorphs (Walls 
1998). Generally, urodelan larvae raised under high-den¬ 
sity conditions took longer to metamorphose and did so 
at a smaller size. It was shown that when A. opacum larvae 
were raised under high-density conditions, which results 
in smaller post-metamorphs, growth and first reproduc¬ 
tion were delayed by three years (Taylor & Scott 1997). 

In the present study, it was shown that although neither 
competition nor sibling-predation (i.e., cannibalism) were 
possible (as potential cannibal larvae were removed once 
identified as such), and food was constantly available, den- 
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Table 3. Density (cm 2 per larva) data from the literature for larvae of different salamander species raised on an unlimited food supply 
(ad libitum). Total numbers of larvae are given in parentheses. In order of appearance data were taken from: Degani (1993), Reques 
& Tejedo (1996), Cohen et al. (2006), Collins & Cheek (1983), Petranka (1984), Wakahara (1995). 


Species 

low density 

medium density 

high density 

food type 

Salamandra salamandra 

75 (20) 

43.7-75 (20) 

7.5 (20) 

aquatic invertebrates 

Salamandra salamandra 

300 (2) 



pond zooplankton 

Salamandra salamandra 

28.7 (4-5) 


14.3 (9-10) 

Tubifex and minced liver 

Ambystoma tigrinum nebulosum 

502.7 (3) 

215.4 (7) 


Artemia nauplii 

Ambystoma texanum 

342 (1) 

85.5 (4) 

42.7 (8) 

aquatic invertebrates 

Hynobius retardatus 

16.5-33 (5) 

13.2-16.5 (20) 

8.25-11 (40) 

urodelan larvae 


sity affected both age and size at metamorphosis. This sup¬ 
ports hypothesis 2 of Petranka & Sih (1986), which states 
that larval growth was density-dependent because of ‘in¬ 
terference competition, a hypothesis that did not find sup¬ 
port in their own study on A. texanum larvae. A similar 


situation existed in the experiment conducted by Eitam et 
al. (2005). The true meaning of this term needs yet to be 
clarified. It is the behavioural aspects of density effects that 
are largely unknown. What could be the exact nature of 
such interference? This term could comprise a number of 
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Figure 7. Minimum and maximum (A) age, (B) mass, and 
(C) length at metamorphosis in the two cohorts. 
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Figure 8. Range between minimal and maximal (A) age, (B) mass, 
and (C) length at metamorphosis in the two cohorts. 
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Table 4. Differences between low and high density (A%) in number, age, and mass at metamorphosis when larvae of different sala¬ 
mander species at low and high densities were raised and fed ad libitum. 


species 

number of larvae (A%) 

age (A%) 

mass (A%) 

source 

Ambystoma opacum 

981 (81.7 %) 

2.6 % 

47.3 % 

Scott (1990) 

Ambystoma opacum 


38.8 % 

49.2 % 

Scott (1994) 

Salamandra infraimmaculata 

5 (50 %) 

4.4 % 

4.9 % 

this study 


factors: (1) Since food (especially live food) is generally not 
evenly distributed, it is possible that some larvae are quick¬ 
er than others in locating it. These will have an advantage 
as they may keep others from feeding. (2) It is possible that 
a sensory response (optical, olfactory or tactile) is involved 
that is sufficiently deterrent to keep some larvae from feed¬ 
ing regularly. Thus, it is theoretically possible that in spite 
of the fact that food was constantly available, some larvae 
fed less or less frequently than others. 

In A. opacum, density affected both larval period (i.e., 
age at metamorphosis), which varied from 2.6 to 38.8%, 
and mass at metamorphosis, which differed by about 48% 
when raised at low and high densities (Scott 1990,1994). 





Figure 9. Average of (A) age, (B) mass, and (C) length at meta¬ 
morphosis in cohorts I and II. See Appendix 1. 


In S. infraimmaculata , the differences between both age 
and mass when raised at low and high densities were 4.4 
and 4.9% respectively (Table 4). 

Differential growth rates of larvae have previously been 
demonstrated in this species (Cohen et al. 2006). These 
were due to different growth modes identified within a 
single half-sibling larval cohort and resulted in variable 
growth rates. Consequently, the variability in the age of 
larvae at metamorphosis could range from two to four 
months. Such wide variability could have catastrophic ef¬ 
fects on the timing of metamorphosis and subsequent ju¬ 
venile dispersal. In this harsh Mediterranean climate, it 
makes all the difference whether larvae metamorphose 
during winter rather than spring, which can well be too hot 
and dry for the post-metamorphs to survive. 

What could be the cause of the definite effect of density 
(larvae raised in the company of other larvae) under seem¬ 
ingly food resource-independent conditions that was dem¬ 
onstrated here? On the face of it, food could not possibly 
be the cause of this effect since it was always freely avail¬ 
able, but perhaps the presence of other larvae already de¬ 
ters some from feeding ad libitum? This point is of great in¬ 
terest and definitely needs further investigations. The evo¬ 
lutionary significance of these findings could be in show¬ 
ing the way how larvae are capable of spreading out the 
timing of metamorphosis (i.e., age at metamorphosis), and 
thereby spreading out post-metamorphic growth, which 
would result in spreading out the age at maturation. This 
would mean that the aquatic larval population were ca¬ 
pable of regulating the size of the adult terrestrial popula¬ 
tion. Keeping the larval population low in small ponds un¬ 
der unpredictable hydroperiod conditions would improve 
chances of survival in species living under unpredictable 
ambient conditions (Warburg 2007b). 
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Appendix 1 

Results of t-tests of data on cohort I versus those on cohort II. 


Age 



single 

low density 

high density 

T value 

0.330005 

0.231763 

0 

P value 

0.745007 

0.819489 

1 

degree of freedom 

19 

17 

30 

Mass 


single 

low density 

high density 

T value 

0.981452 

1.650592 

1.889448 

P value 

0.343027 

0.111849 

0.068216 

degree of freedom 

14 

24 

31 

Length 


single 

low density 

high density 

T value 

-0.25577 

1.543785 

-0.87974 

P value 

0.801845 

0.135725 

0.385772 

degree of freedom 

14 

24 

31 


Appendix 2 

Results of regression analysis of age under different density 
regimes in relation to Figures 1 and 2. 


Cohort I 



R 2 

P 

degree of freedom 

single versus age 

0.9313 

0.531 

3 

low density versus age 

0.7767 

0.2889 

4 

high density versus age 

0.757 

0.806 

3 

Cohort II 


R 2 

P 

degree of freedom 

single versus age 

0.8401 

0.0598 

9 

low density versus age 

0.7281 

0.7432 

13 

high density versus age 

0.4973 

0.8015 

14 


Appendix 3 

Results of regression analysis of age under different density 
regimes in relation to Figures 3 and 4. 


Cohort I: single 



age 

mass 

length 

R 2 

0.9941 

0.9981 

0.9981 

P 

0.5193 

0.1829 

0.1829 

degree of freedom 

6 

6 

6 


Cohort I: low density 



age 

mass 

length 


Age 

Mass 

Length 

R 2 

0.9779 

0.9989 

0.9996 

P 

0.5879 

0.0001 

< 0.0001 

degree of freedom 

20 

20 

20 

Cohort I: high density 


age 

mass 

length 


Age 

Mass 

Length 

R 2 

0.9779 

0.9982 

0.9989 

P 

0.6382 

0.0261 

< 0.0001 

degree of freedom 

17 

17 

17 

Cohort II: single 


age 

mass 

length 

R 2 

0.9963 

0.095 

2 

P 

0.9993 

< 0.0001 

2 

degree of freedom 

0.9986 

0.0198 

2 

Cohort II: low density 


age 

mass 

length 

R 2 

0.983 

0.9984 

1 

P 

< 0.0001 

< 0.0001 

0 

degree of freedom 

16 

16 

16 

Cohort II: high density 


age 

mass 

length 

R 2 

0.9974 

0.9984 

0.9982 

P 

0.1882 

0.0703 

0.0085 

degree of freedom 

10 

10 

10 


Appendix 4 

Results of regression analysis of age under different density 
regimes in relation to Figures 5 and 6. 


Cohort I 



R 2 

P 

degree of freedom 

age versus mass 

0.9941 

0.5193 

6 

age versus length 

0.9981 

0.1829 

6 

mass versus length 

1 

0.06 

6 

Cohort II 


R 2 

P 

degree of freedom 

age versus mass 

0.9974 

0.1882 

10 

age versus length 

0.9984 

0.0703 

10 

mass versus length 

1 

0 

10 
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